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Abstract: Anti-leukemia screening of previously prepared isothiouronium and quaternary 
salts was performed, and some salts exhibited promising activity as anticancer agents. 
Quantum chemical calculations were utilized to explore the electronic structure and 
stability of these compounds. Computational studies have been carried out at the PM3 
semiempirical molecular orbitals level, to establish the HOMO-LUMO, IP and ESP 
mapping of these compounds. The ADMET properties were also studied to gain a clear 
view of the potential oral bioavailability of these compounds. The surface properties 
calculated included critical micelle concentration (CMC), maximum surface excess (Fmax), 
minimum surface area (Amin), free energy of micellization (AG^mic) and adsorption (AG^ads)- 
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1. Introduction 

Cancer is presently considered as one of the primary worldwide health problems, responsible for 
about 15% of deaths and 25% in developed countries. Antitumor chemotherapy is thus nowaday s a 
very attractive research target [1,2]. Thioureas, isothiouronium compounds and their derivatives 
constitute an important class of compounds which exhibit a wide range of antibacterial, fungicidal. 
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herbicidal, antiviral, and plant growth regulatory activities, and play important roles in many chemical 
and biological processes [3-7], but also act as potential anticancer and anti-HIV drugs [8-11]. The 
activity of isothiouronium groups may be due to the enhanced acidity of the NH moieties, thereby 
functioning as a better binder than the thiourea group [12-17]. Several synthetic isothiouronium 
compounds were synthesized to obtain more potent and less toxic therapeutic agents [18]. The studies 
report herein studies aimed to examine the anti-cancer and ADMET properties of some previous 
synthesized isothiouronium salts 1^ (Figure 1) [19]. 

Figure 1. The chemical structure of the previously prepared isothiouronium salts [19]. 




2. Result and Discussion 

2.7. Surface Parameter Results for the Prepared Isothornium Salts Surfactants 
2.1.1. CMC of the Prepared Surfactants 

Surfactants form aggregates of molecules or ions called micelles, when the concentration of the 
surfactant solute in the bulk of the solution exceeds a limiting value, the so called critical micelle 
concentration (CMC), which is a fundamental characteristic of each solute-solvent system. If the 
properties of a surfactant solution are plotted as a function of the concentration of the surfactant, the 
properties usually vary linearly with the concentration, up to critical micelle concentration, at which 
point there is a break in the curve as shown in Figure 2. 

The results in Figure 1 and Table 1 showed that all tested compounds have very low CMC values 
and a notable decrease of these values was observed in going from 1 to 4. That fact could be explained 
from the unique property of the quaternary ammonium salts in water, that is the salts retain their unity 
in their solutions, which increases their volume in the aqueous media and then repulsion occurs 
between the hydrophobic chains and water molecules. The lower CMC values of 2 compared with 1 
and 4 compared with 3 is due owing to the same reason, whereby an increased hydrophobic part 
increases this repulsion and accelerates the adsorption of molecules at the interface under very low 
concentration that form micelles immediately. 



Pharmaceuticals 2013, 6 



636 



Figure 2. Variation in surface tension of surfactants 1-4 vs. concentration at 25 °C 
concentration, expressed as mol/L. 
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Table 1. The critical micelle concentration (CMC) and surface parameters of the tested 
compounds (surfactants). 



Comp. 


CMC 


YCMC 

(mN/m) 


IIcMC 

(mN/m) 


Pc20 

(Mole/L) 


Tmax X 10 

(Mole/cm^) 


A min 

(nm^) 


^ Gads 


A Gmic 


1 


1.3 


51 


21 


3.0 


6.9 


2.3 


-36 


-33 


2 


1.05 


49 


23 


3.2 


6.5 


2.5 


-37.5 


-33.9 


3 


0.9 


41 


31 


4.5 


4.9 


3.3 


-41 


-34.7 


4 


0.6 


38 


34 


5.0 


4.5 


3.6 


-44 


-36.4 



2.1.2. Effectiveness (ITcmc) 

The effectiveness (ITcmc) is determined by the difference between surface tension values at the 
CMC "ycmc" and the surface tension measured for pure water at the appropriate temperature (yo). The 
most efficient surfactant one is that gives the greatest lowering in surface tension for a critical micelle 
concentration (CMC). Compound 4 was found to be the most efficient one in Table 1 because it 
achieved the maximum reduction of the surface tension at CMC (Figure 2). 

2.1.3. Efficiency (PC20) 

The efficiency (PC20) was determined by the concentration (mol/L) capable of suppressing the 
surface tension by 20 dyne/cm. Values of efficiency of the prepared surfactants are shown in Table 1, 
the efficiency increased from compound 1 to 4. 

2.1.4. Maximum Surface Excess (Fmax) 

The number of surfactant molecules at the air-water interface at the critical micelle concentration at 
(25 "^C) is expressed by Fmax. It is evident from Table 1 that Fmax decreases from 1 to 4 this due to the 
increase in size of the molecules, which are adsorbed at the interface by its huge size and then form the 
micelles directly at low concentration. 
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2.1.5. Minimum Area per Molecule (Amin) 

The minimum surface area is defined as the area occupied by a surfactant molecule at the air-water 
interface at the equilibrium of the solution. Results given in Table 1 indicate that the decrease of Fmax 
that occurred at the interface was logically accompanied by an increase in Amin values from 1 to 4. 

2.1.6. The Standard Free Energies of Micellization (AG Vic) and Adsorption (AG^ads) 

From Table 1 the values of AG^mic and AG^ads are always negative, indicating the spontaneity of 
these two processes, although there is more increase in the negativity of AG^ads than that of 
micellization, indicating the tendency of the molecules to be adsorbed at the interface. 

2.2. Anti-Cancer Screening 

Tumor growth can be checked by blocking the process of multiplication of tumor cells; this could 
be achieved either by irradiation or chemical agents. The anticancer activity of the tested compounds 
1-4 on the life-span of mice bearing Leukemia L-1210 was evaluated (Table 2). 

Table 2. Anticancer activity of some selected compounds against leukemia L-1210. 



Compds. 


Dose (fig/kg.) 


ILS % * 




400 


0 


1 


200 


100 




100 


100 




400 


0 


2 


200 


91 




100 


99 




400 


0 


3 


200 


105 




100 


120 




400 


61 


4 


200 


107 




100 


114 



* Percentage ILS (increase in life span) values were determined from median survival times of 
treated and saHne control groups. If more than 50% of the animals survived for more than 30 days, 
median survival times and % ILS are indicated as greater than 30 days. 

As shown, the isothiuronium compounds 1 and 2 were toxic at a dose level of 400 |ig/kg, and failed 
to increase the life-span of the treated animals. In addition, the quaternary Entobex derivatives 3 and 4, 
exhibited borderline activity amounting to 20 and 14%, respectively. From the combination between 
surface and antitumor activity we reach the conclusion that compounds 3 and 4 with the highest 
surface activity were also the compounds that exhibited antitumor activity. 

In our research surfactants 3 and 4 affected tumor tissues at very low concentrations, lower than 
their CMC values, which means that there is a strong relation between very small CMC values of these 
compounds and the display of borderline activity at very low concentration, this being due to the fact 
that an increasing concentration of surfactant causes an increase in the adsorption process on cell 
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membranes until the CMC is reached, and beyond this concentration the adsorption is slowly retarded 
and finally stopped due to the formation of micelles which prevent the mobility and suppress 
anticancer activity. 

2,3, Molecular Modeling 

In trying to gain better insight into the molecular structure of compounds 1-4 in Figure 1, conformational 
analysis of the target compounds has been performed using the MMFF94 force-field [20,21] 
(calculations in vacuo, bond dipole option for electrostatics, Polake Ribiere algorithm, RMS gradient 
of 0.01 kcal/A mol) as implemented in the Spartan 08 program [22]. The most stable conformer was 
fully geometrically optimized by PM3 [23] semi-empirical Hamiltonian molecular orbital calculation 
[with Restricted Hartree-Fock (RHF) and RMS gradientof 0.05 kca]. Furthermore, the computed 
molecular parameters, total energy, heat of formation, the lowest occupied molecular orbital (LUMO) 
and the highest occupied molecular orbital (HOMO) energies, potential energies, solvation energies, 
electrostatic energies, ionization potential energies, dipole moment and the thermo-dynamic 
parameters of the studied compounds, such as zero-point vibrational energies, entropy, constant 
volume molar heat capacity, enthalpy and Gibbs free energy, were calculated (Tables 3 and 4). 



Table 3. The optimized calculatedenergies and thermodynamic parameters atthe semi-empirical 
PM3 level for 1-4. 



CPD 


HF 


E 


Es 


ZPE 


S° 


Cv° 


H° 


G° 


1 


-458.26 


-407.84 


-43.58 


1601.57 


462.36 


236.63 


0.155427819 


-0.207933 


2 


-376.38 


-446.96 


-70.57 


628.37 


516.22 


296.04 


0.112537693 


0.0539166 


3 


-1835.03 


-1914.11 


-79.08 


514.11 


606.16 


385.66 


-0.436475908 


-0.505311 


4 


-508.01 


-628.05 


-45.80 


602.05 


857.66 


697.30 


0.452964763 


0.3555698 



HF: The heat of formation energy (kJ/mol); E: Total energy (kJ/mol) (heat of formation + strain energy), 
Es: Solvation energy (kJ/mol), ZPE: zero-point vibrational energies (kJ/mol), S°: Entropy (kJ/mol), 
Cv'^: Constant volume molar heat capacity, H°: Enthalpy (kJ/mol), G°: Gibbs free energy (kJ/mol). 



Table 4. The optimized calculatedenergies atthe semi-empirical PM3 level for 1-4. 



CPD 


HOMO 


LUMO 


HLG 


IP 


ESP 




Z 


1 


-10.21 


-2.31 


7.9 


299.90 


77.90 


6.04 


3.95 


2 


-8.66 


-1.32 


7.34 


51.61 


-285.32 


5.74 


3.67 


3 


-8.03 


-1.95 


6.08 


43.46 


-448.58 


9.84 


3.04 


4 


-8.77 


-1.51 


7.26 


42.94 


-279.41 


7.97 


3.63 



HOMO: Highest Occupied Molecular Orbital (eV); LUMO: Lowest unoccupied Molecular Orbital (eV); 
HLG: difference between HOMO and LUMO energylevels; IP: ionization potential energy (kJ/mol); 
ESP: electrostatic potential energy (kJ/mol), /: Mulliken electronegativity (eV), ji: dipole moment (Debye). 

2.3.1. Energy Features 

The calculated molecular parameters have been used to investigate the most stable conformer of the 
compounds 1-4 (Table 3), and the optimized geometry structures are represented in Figure 3. The 
lower calculated energy of the most active compounds 3 and 4 (-1914.11 and -628.05 kJ/mol) 
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compared with the inactiveisothiouronium salts 1 and 2 (-407.84 and -446.96 kJ/mol, respectively), 
suggests an increase in spatial stability for 3 and 4 compared to 1 and 2. Furthermore, compound 3 has 
the lowest calculated solvation energy (-79.08 kJ/mol) which may explain the increasing activity due 
to increased solubility. Taking into consideration the zero-point energy at 298.15 K temperature, the 
inter-conversion energy is the largest for 1 and 2. From comparison of the thermodynamic parameters 
(Table 4), the most active compounds 3 and 4 are the more stable structures, which agrees with the 
PM3 calculations. 

2.3.2. Optimized Geometry 

The single crystal x-ray structures of the isothiouranium compounds 1-4 are not available till now, 
so the geometrical optimization parameters (bond lengths and angles) at thePM3 semi-empirical 
quantum mechanical level (Figure 3, Figure SI in Supplementary Information) were calculated 
(Tables S1-S5 in Supplementary Information). 

Figure 3. Ball and stick rendering for the most stable conformer of the isothiouronium 
salts 1-4 as calculated by PM3 semi-empirical molecular orbital calculations. 



1 2 




In 1 and 2, the N1-C2-N3 and N9-C8-N10 bond angles are ~1 19°, and the two carbamimidothioate 
moieties are stabilized by being arranged in a coplanar position with each other, and the distances 
between C2-N3, Nl-C2and C2-S4 are nearly the same as N10-C8, C8-N9 and C8-S7(^1.3, ^1.409 and 
-1.8 A) respectively, and the increased bond length of the imino moiety could be explained by the 
presence of tt electron conjugation of the phenyl ring, and the sp hybridization N-imino orbitals that 
has more s-character, with the most electron density closer to the nuclei compared with the electron 
density distribution at sp^ hybridized amino nitrogen, which results in a lengthening of the C2-N1 and 
C8-N9 imino bond. The compounds 3 and 4 are arranged in space in the same manner, except the 
methyl moiety is arranged perpendicular with the phenanthrolinedione ring in 3, and in 4, the propionic 
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moiety is arranged coplanar with the phenanthroUnedione ring. The bond lengths in both compounds 3 
and 4 for C1-C6, N14-C8, N14-C13 and C9-C10 are about ^1.33, 1.4, 1.31 and l.SOOA, respectively. 

2.3.3. HOMO-LUMO and Dipole Moments Analysis 

The frontier molecular orbitals (highest occupied molecular orbital, HOMO, and lowest unoccupied 
molecular orbital, LUMO) are the most important orbitals in a molecule. These orbitals determine the 
way the molecule interacts with other species. The frontier orbital gap (as calculated by simple Hiickel 
Molecular Orbital theory, SHMO) helps to characterize the chemical reactivity and kinetic stability of 
amolecule [24-26]. Furthermore, the HOMO and LUMO of a molecule play important roles in 
intermolecular interactions [27], through the interaction between the HOMO of the drug with the 
LUMO of the receptor and vice versa. The interactions are stabilized inversely with the energy gap 
between the interacting orbitals. Increasing HOMO energy and decreasing LUMO energy in the drug 
molecule lead to enhancement of stabilizing interactions, and hence, binding with the receptor [27]. 
HOMOs and LUMOs of the studied systems in the So states are shown in Figure 4, which suggests the 
delocalization and localization of molecular orbitals. 

Figure 4. The calculated HOMO, LUMO and HOMO-LUMO gap for 1-4 as calculated by 
the PM3 semiempirical method. 



The HOMO/LUMO and gap values for compounds 1-4 are given in Table 4. The results showed 
that 3 and 4 have the lowest energy gaps (7.26 and 6.08 eV), and explains their higher affinity towards 
cancer cells. The total dipole moment is also considered an important physical quantity, as it reflects 
the interaction ability of the molecules with the surrounding environment. Compounds 3 and 4 have 
higher dipole moment values (-9.84 and 7.97 Debye) than 1 and 2 (-6.04 and 5.74 Debye), 
respectively, which increases their ability to interact with the surroundings. Therefore, molecules 3 and 
4 are promising structures for interaction with the receptors of cancer cells. 
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2.3.4. Flexible Alignment 

To understanding the similarity between the three-dimensional structures of active compounds 3 
and 4, flexible alignment using MOE/MMFF94 was employed (Figure 5) [28]. 

Figure 5. (a) Flexible alignments of the most active compounds 3 (in green), 4 and (in 
red), (b) flexible alignment of the highly active compound 3 (in red) and lowest active 
compounds 1 (in yellow) and 2 (in blue). 




Two hundred conformers of each compound were generated and minimized with a distance-dependant 
dielectric model. A low energy set of 100 was selected for further analysis. Conformations of compound 
3 were generated using distance geometry and optimized with MMFF94. Thirteen low energy, 
maximally dissimilar structures were selected for comparison to the other compounds. After assigning 
MMFF94 charges to all molecules, flexible alignment was ranked by overlays of compounds 3 and 4 
based on electrostatic, steric field, hydrophobic areas overlap, hydrogen bond acceptors and donors 
overlap. From the highest scoring superposition, the limited set of conformers was used in the analysis 
of molecules with high flexibility capable to achieving complete atom to atom superposition. 

A common feature of the MOE-generated alignments was that the two structures showed matched 
phenanthrolinedione rings, and slight differences in the distance the between methyl and acetyl 
moieties (0.28A, Figure 4a). In the same way, the most active anti-cancer compound 3 and the least 
active compounds 1 and 2 were subjected to flexible alignment. Analysis of the weakly active 
molecules (Figure 4b), is an important way to gain a clear vision of the essential features for a given 
activity. It is clear that compounds 1 and 2 were flexibly aligned in a different manner when compared 
with the active compound 3. Common feature of the MOE-generated alignments showed, the super 
position of the dicarbamimidothioate fragments in 1 and 2, and a deviation distance of about 2.20A 
compared with the highly active compound 3. These features explaining why these compounds were 
the least active as anti-cancer agents and shows the importance of both the alkyl moiety and the 
heteroaryl ring attached to it (Figure 5). 

2.3.5. Electrostatic (ESP) and Ionization Potential (IP) Map 



In an attempt to understand the lowest and the highest anti-cancer activity of compounds 1-4, 
electrostatic and ionization potential mapping was performed for the lowest energy conformers, to 
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examine the match and mismatch in electronic, and electrostatic binding characteristics for the 
molecules surface and conformational properties (Figure 6). 

Figure 6. The ESP and IP surfaces for 1-4 by PM3 semi-empirical calculations. 



ESP 



IP 



The mapped electrostatic potentials (ESP), represent a balance between the repulsive interactions of 
the positively-charged nuclei and attractive interactions involving the negatively-charged electrons. 
The colors toward red depict negative potential (high electron density area, representing a strong 
attraction between the proton and the points on the molecular surface), while colors toward blue depict 
positive potential, and colors in between (orange, yellow, green) depict intermediate values of the 
potential. In addition, the local ionization potential (IP) mapping was performed, which indicates the 
proportional electron ionization facility around the molecule [29,30] and measure the sensitivity of a 
molecular zone toward electrophilic attack (reactivity). The default program value used [Fixed Value: 
(0.032e); Medium Resolution], for the IP surface mapping is shown in Figure 6. Comparison of the 
ESP and IP mappings of the lowest activity compounds 1 and 2 showed a common feature of 
increasing distribution of positive charge (in blue), and increase the polar area. On the contrary, the 
high activity compounds 3 and 4 showed an increased distribution of negative charge (in red) on the 
surface of the compounds. From the above features, the highest activity may be due to the increasing 
negative charge located on the compounds, representing hydrogen bond acceptors (in red). It is clear 
that the related charge distribution and electrostatic ionization mappings suggest a similar interaction 
and orientations of the molecules with a potential protein-binding site. 

2,4. ADMET Factor Profiling 

Oral bioavailability is considered to play an important role in the development of bioactive 
molecules as therapeutic agents. Many potential therapeutic agents fail to reach the clinic because of 
their unfavorable absorption, distribution, metabolism, elimination and toxic (ADMET) factors. 
Therefore, a computational study for prediction of ADMET properties was performed for compounds 
1-4, by the determination of topological polar surface area (TPSA), a calculated percent absorption 
(%ABS) which was estimated by the Zhao et al. equation [31], and the "rule of five", formulated by 
Lipinski [32], which establish whether a chemical compound can be expected to be an orally active 
drug in humans, which occurs if it has no more than one violation of the following rules: (i) ClogP 
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(partition coefficient between water and octanol) < 5; (ii) number of hydrogen bond donors sites < 5; 
(iii) number of hydrogen bond acceptor sites < 10; (iv), molecular weight < 500 and molar refractivity 
should be between 40-130. In addition, the total polar surface area (TPS A) is another key property 
linked to drug bioavailability, a the passively absorbed molecules with TPS A > 140 have low oral 
bioavailability [33]. All calculated descriptors were obtained using the MOE package [32], and the 
results are listed in Table 5. 

Table 5. Pharmacokinetic parameters important for good oral bioavailability of compounds 1-4. 



CPD 


Mwt 


TPSA 


%ABS 


V 


LogP 


HBD 


HBA 


Lip-V 


1 


408.344 


187.58 


44.2849 


255.875 


1.58 


4 


8 


0 


2 


642.76 


249.88 


22.7914 


658.961 


2.27 


4 


14 


2 


3 


352.13 


50.88 


91.4464 


265.068 


2.83 


0 


4 


0 


4 


363.16 


88.50 


78.4675 


307.261 


2.035 


1 


6 


0 



TPSA: Polar surface area; %ABS: absorption percentage; Log P: Calculated lipophilicity; LogS: Solubility 
parameter; HBA: Number of hydrogen bond acceptors; HBD: Number of hydrogen bond donors; Lip-V: 
Number of violations of Lipinski's rule of five; V: Volume (A^). 

Our results revealed that the CLogP (the lipophilicity factor) [34] w^as less than 5.0, the molecular 
v^eight was less than 500, except for compound 2 (642.76), hydrogen bond acceptor < 10 except for 2 
(14) and hydrogen bond donors (-1^), which fulfill Lipinski's rule. Also, the percent absorption of 
compounds 3 and 4 show the highest % absorption among all compounds 1-4. From these data, we 
can suggest that 3 and 4 can be used as good orally absorbed anti-cancer compounds with less toxicity, 
but compound 2 is highly toxic and has the lowest oral absorption from all investigated compounds, 
which matches the empirical data. 

2.5. Pharmacophore Prediction 

The aim of this approach is generation and prediction of a pharmacophore model (hypothesis) based 
on the most active compound. The following steps are employed for 3D pharmacophore applications: 
(i) Compound 3 is used as reference for other conformations of each compound (Figure 6); (ii) 
Assignment of the pharmacophoric features; (iii) Application of the method for undertaking 
conformational searching of databases for new structures matching the generated pharmacophoric 
features. These steps allow (hypothetical) pharmacophoric building of anticancer compounds. MOE is 
used for pharmacophore building [28]. The calculated conformational models were performed with a 
15 Kcal energy cut off (minimization convergence criteria during conformational analysis: energy 
convergence = 0.01 kcal/mol, gradient convergence = 0.01 kcal/mol). The number of conformers 
generated for each substrate was limited to a maximum of 500. All molecules with their associated 
conformations were regrouped depending on the biological data. Hypothesis generation was carried 
out with the low energy conformers of the molecules. The calculation and analysis were performed 
after assignment of possible pharmacophore elements for each analogue using the MOE program, then 
superposition of the molecules including the assigned elements was attempted (Figure 7), to generate 
pharmacophore maps, and many runs of the calculations were repeated. 
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Figure 7. (a) Showsthe best predicted pharmacophore features and geometries which are 
required for anti-cancer activity; (b) showsthe most active compound 3 mapped to the 
pharmacophore model; (c) showsthe least active compound 2 mapped to the pharmacophore 
model. Pharmacophore features are color coded: yellow for catalytic aromatics, green for 
hydrophobic, red for acceptor aromatics and blue a hydrogen bond acceptor feature. 




For each run a distinct number of specified pharmacophore elements were generated. All 
appropriate models showed that the acceptor atoms of the CO fragments (as acceptor) and the methyl 
fragment (as hydrophobic unit) were well superimposed within the set distance tolerance, which 
confirms the important role of the hydrophobic moieties for anticancer activity. The initial pharmacophoric 
query was carried out, and the introduction of four features, as summarized in Table 6, is illustrated 
graphically in Figure 7. 

Table 6. Pharmacophoric and structure features of the most active compounds. 



Pharmacophoric features 


Structure features 


Fl: Arc/Cat 


Pyridinyl 


F2: Aro/Acc 


Alkyl pyridinyl 


F3: Ace 


CO 


F4: Hydro 


Alkyl 



From Figure 7a, models for 3 and its analog 4 possess pharmacophore elements: the heteroaryl ring 
(as a hydrophobic moiety) attached to an alkyl group which is in plane with the other heteroaryl ring, 
two heteroaryl rings attached with a CO group (as a hydrogen bond acceptor region. Figure 7b). In 
contrast, the model for compound 1 and its analog 2 (Figure 7c) showed mismatch with the mapping of 
the pharmacophore model (Figure 6a), since the alkyl attached with the heteroaryl part of compound 3 
plays an important role in activity. According to the pharmacophore generated by MOE [28] the 
minimal structural requirements for anti-cancer activity consist of an alkyl fragment (hydrophobic 
region) attached to a heteroaryl ring (aryl catalysis region), carbonyl fragment (H-bonding acceptor 
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region) and another hetroaryl ring (aryl acceptor region). This pharmacophoric assumption was 
consistent with biological data. 

3. Experimental 

3.1. Surface Active Properties Evaluation Methods 

3.1.1. Surface Tension Measurements 

Surface tension measurement of the prepared surfactants was made at (25 "^C) with Du 
Nouytensiometer (Kruss type 8451) using distilled water solution of 0.1% weight concentration [35]. 
The surface tension of the used distilled water was 73 mN/m, surfactant solutions were aged for 1/2 h 
before any measurements were made. Three readings were made on each sample to determine any 
change with time and to obtain an average value [36]. 

3.1.2. Efficiency (PC20) 

The efficiency (PC20) was determined as the concentration (mol/L) capable of suppressing the 
surface tension by 20 dyne/cm [37]. The efficiencies have been determined by extrapolating from 
y = 52 to the linear portion before CMC of the y versus-Logc plot [38], at 25 °C. 

3.1.3. Effectiveness (Ilcmc) 

The surface tension (ycmc) values at CMC were used to calculate values of surface pressure 
(effectiveness) from the following expression: 

ncmc ~ yo ~ ycmc 

where yo is the surface tension measured for the pure water at the appropriate temperature and ycmc is 
the surface tension at the CMC. The effectiveness of adsorption is an important factor to determine 
properties of a surfactant such as foaming, wetting and emulsification, since tightly packed coherent 
interfacial films have very different interfacial properties than loosely packed, non-coherent films [39]. 

3.1.4. Determination of Critical Micelle Concentration 

CMC of the prepared surfactants were determined by the surface tension method [40]. In this 
method, values of the surface tension obtained for various concentrations of aqueous solutions of the 
prepared surfactants were plotted v^*. the corresponding concentrations. 

3.1.5. Maximum Surface Excess Fmax 

The surface excess concentration is defined as the surface concentration at surface saturation; the 
maximum surface excess Fmax is a useful measure of the effectiveness of adsorption of the surfactant at 
the water-air interface, since it is the maximum value at which adsorption can attain: 

r 1 r -I 

^^^^ 2.303RT^51ogc^T 
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where R = 8.314 J mol"^ K~\ T is absolute temperature, (Sy/S logc) is the slope of the surface tension 
v^*. log concentration plot at 25 °C [41]. A substance which lowers the surface tension is thus present in 
excess at or near the surface, i.e., when the surface tension decreases with increasing activity of the 
surfactant, F is positive. 

3.1.6. Minimum Surface Area (Amm) 

Amin is the minimum area per molecule of the prepared compounds at the interface and was 
calculated from the following equation: 



where N is Avogadro's number and Fmax is the maximum surface excess. 

3.1.7. The Standard Free Energies of Micellization AGVic and Adsorption AG^ads 

Understanding the process of micellization and adsorption are important for explanation of the 
effects of structural and environmental factors on the value of the CMC and for predicting itseffects on 
new structural and environmental variations. Standard free energy of micellization AG^mic, and 
adsorption AG^ads, have played an important role in such understanding.The standard free energy of 
micellization and adsorption are given by: 

AGVic = RTln CMC 
AG^ds = AGVic - 6.023 x IO-^Hcmc A^^in 

3.2. Anticancer Screening against Leukemia L-1210 

Group of six animals were inoculated interperitoneally with 10^ cells. After 24 h the tested 
compounds were injected at three dose levels selected on the basis of their toxicity so that the higher 
dose is the maximum treatment dose. The animals were kept for one month or till the death of the last 
animal. The ILS were calculated for all animals died from tumors, and T/C values were calculated. 

3.3. Molecular Modeling 

3.3.1. Computational Details 

All electronic structure calculations were performed using the Spartan 08 program [22]. Geometry 
optimizations have been achieved using the PM3 method. The structural parameters, such as the dipole 
moment of the molecules the energy of the highest occupied molecular orbital EHOMO and the 
lowest unoccupied molecular orbital ELUMO were obtained. The thermodynamic parameters, 
including total energy(TE), enthalpy (H), entropy (S), free energy (G), zero-point vibrational energy 
(ZPE), the absorption spectra and the vibration frequency calculation For each stationary point 
was performed at the same levels, to characterize its nature as minima or transition states and to 
correct energies. 
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3.3.2. Conformational Search and Flexible Alignment 

Conformational analysis and flexible alignment of tested compounds were carried out with the 
MOE [28] software of Chemical Computing Group Inc. (Montreal, Canada), on a Core 2 duo 3.00 GHz 
workstation. The molecules were built using the Builder module of MOE. The geometry optimized 
using the MMFF94 force-field followed by a flexible alignment with systematic conformational 
search. Lowest energy aligned conformation(s) were identified through the analysis module of DSV. 

4. Conclusions 

Anticancer screening of some previously prepared novel isothiouronium and quaternary salts was 
performed, and compounds 3 and 4 showed promising activity as anti-leukemia agents. Computational 
studies have been performed at the PM3 semiempirical molecular orbitals level, to establish the 
HOMO-LUMO, IP and ESP mapping of these compounds, to understand the lowest and highest 
anticancer activity of the compounds, and show the importance of the alkyl group in these compounds. 
Also, the ADMET properties for compounds 3 and 4 showed good oral absorption and suggested they 
could be used as anti-cancer compounds with less toxicity. 
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